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Abstract 

The  start  up  behaviour  of  PEM  fuel  cells  below  0  °C  is  one  of  the  most  challenging  tasks  to  be  solved  before  commercialisation.  The  automotive 
industry  started  to  develop  solutions  to  reduce  the  start  up  time  of  fuel  cell  systems  in  the  middle  of  the  nineties.  The  strategies  varied  from 
catalytic  combustion  of  hydrogen  on  the  electrode  catalyst  to  fuel  starvation  or  external  stack  heating  via  cooling  loops  to  increase  the  stack 
temperature. 

Beside  the  automotive  sector  the  cold  start  ability  is  as  well  important  for  portable  PEMFC  applications  for  outdoor  use.  But  here  the  cold  start 
issue  is  even  more  complicated,  as  the  fuel  cell  system  should  be  operated  as  passive  as  possible. 

Below  0  °C  freezing  of  water  inside  the  PEMFC  could  form  ice  layers  in  the  electrode  and  in  the  gas  diffusion  layer.  Therefore  the  cell  reaction  is 
limited  or  even  inhibited.  Product  water  during  the  start  up  builds  additional  barriers  and  leads  to  a  strong  decay  of  the  output  power  at  isothermal 
operating  conditions. 

In  order  to  find  out  which  operational  and  hardware  parameters  affect  this  decay,  potentiostatic  experiments  on  single  cells  were  performed 
at  isothermal  conditions.  These  experiments  comprise  investigations  of  the  influence  of  membrane  thickness  and  different  GDL  types  as  well  as 
the  effect  of  gas  flow  rates  and  humidification  levels  of  the  membrane.  As  pre  stage  to  physical  based  models,  empirical  based  prediction  models 
are  used  to  gain  a  better  understanding  of  the  main  influence  parameters  during  cold  start.  The  results  are  analysed  using  the  statistical  software 
Cornerstone  4.0. 

The  experience  of  single  cell  investigations  are  compared  to  start  up  behaviour  of  portable  fuel  cell  stacks  which  are  operated  in  a  climate 
chamber  at  different  ambient  temperatures  below  0  °C.  Additional  flow  sharing  problems  in  the  fuel  cell  stack  could  be  seen  during  cold  start  up 
experiments. 
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1.  Introduction 

Fuel  cells  are  known  as  one  of  the  most  promising  energy 
conversion  systems  for  the  future.  Their  principle  of  direct  con¬ 
version  of  chemical  energy  into  electrical  energy  and  heat  let 
them  fit  to  various  applications.  Especially  the  automotive  indus¬ 
try  pushed  the  PEMFC  development  in  the  recent  years  strongly. 
Beside  very  high  dynamic  demands  in  cars  fuel  cell  systems  must 
also  be  able  to  start  at  temperatures  below  0  °C. 

Operational  demands  in  other  fuel  cell  applications  are  nor¬ 
mally  much  lower.  For  portable  applications  with  a  fuel  cell 
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system  power  below  100W  a  lot  of  promising  systems  were 
presented  which  show  good  performance  data  at  ambient  con¬ 
ditions  [1-4]. 

But  especially  for  the  advanced  hydrogen/air  PEMFC  tech¬ 
nology  new  challenges  come  up.  The  ability  to  deliver  power 
to  outdoor  applications  like  independent  power  sources  for  out¬ 
door  environment  or  traffic  gauging  stations  demands  cold  start 
ability  of  a  fuel  cell  system. 

Compared  to  the  automotive  fuel  cell  systems  portable 
PEMFC  systems  should  be  more  simple.  For  cold  start  no  periph¬ 
eral  system  components  should  be  used  with  respect  to  the 
overall  system  efficiency.  The  aim  of  our  investigations  was  to 
find  out  the  main  influencing  parameters  on  the  fuel  cell  cold  start 
ability  covering  operational  parameters  and  hardware  parame¬ 
ters  in  single  fuel  cells  as  well  as  in  portable  fuel  cell  stacks. 
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2.  Experimental  set-up 

2.7.  Single  cell  test  bench 

A  test  bench  was  set-up  which  enables  single  fuel  cell 
operation  at  nearly  isothermal  conditions.  The  temperature 
(between  —16  and  70  °C)  can  be  controlled  by  a  cryostat  bath 
(Lauda)  which  is  connected  to  the  fuel  cell  cooling  plates. 
The  gas  humidity  is  adjustable  up  to  100%  using  conven¬ 
tional  bubble  humidifiers.  The  gas  inlet  and  outlet  pressures, 
dew  points,  and  temperatures  are  monitored  as  well  as  the 
cell  temperature  of  the  anode  and  cathode  graphite  flowfield 
plate.  The  whole  test  bench  is  controlled  automatically  by 
Lab  View®.  All  signals  are  recorded  by  a  data  logger  (Agilent). 
For  gas  regulation  digital  mass  flow  controllers  (Bronkhorst)  are 
used.  The  electronic  load  (Zentro  EL)  has  an  operating  range 
up  to  50  A. 

The  single  test  cell  has  an  active  area  of  46  cm2.  The  flow- 
field  has  a  double  serpentine  structure  with  a  channel  width  of 
1.2  mm,  depth  0.8  mm  and  a  ligament  width  of  1.2  mm.  The 
fuel  cell  assembly,  comprising  the  catalyst  coated  membrane, 
gas  diffusion  layer  (GDL),  graphite  compound  flowfield  plates, 
titanium  end  plates,  and  plastic  cooling  plates  are  placed  in  an 
isolation  housing  made  of  Styrodur®  C  to  minimize  the  impact 
of  heat  exchange  with  the  ambience. 

2.2.  Climatic  chamber  test  bench 

For  realistic  fuel  cell  system  testing  a  test  bench  with  a  cli¬ 
matic  chamber  (Votsch  VC  4034)  was  set-up.  Fuel  cell  systems 
and  fuel  stack  operational  behaviour  can  be  investigated  in  a 
range  between  —40  and  100  °C  and  10%  and  98%  r.h.  All  main 
components  are  controlled  by  Lab  View®  and  all  signals  are  mon¬ 
itored  via  data  logger  and  PCI  measuring  PC  card.  Gases  can  be 
supplied  humidified  or  dry  to  the  fuel  cell.  The  electronic  load 
(Zentro  EL)  can  operate  up  to  125  A. 

The  fuel  cell  stacks  consist  of  six  cells  with  an  active  area 
of  30  cm2  of  each  cell.  Bipolar  plates  were  made  of  graphite 
compound  with  milled  double  serpentine  flowfields.  The  current 
collectors  are  made  of  copper  and  the  endplates  of  aluminium. 

The  fuel  cell  stack  was  placed  in  a  plastic  housing,  to  pre¬ 
vent  convective  cooling  by  the  climate  chamber  fan.  The  air  is 
supplied  by  a  pump  (Nitto)  which  delivers  cold  air  from  the  inte¬ 
rior  of  the  climate  chamber.  Dry  hydrogen  is  supplied  by  the  test 
bench  and  the  gas  temperature  is  reduced  by  cooling  coils  inside 
the  climate  chamber  interior. 

3.  Experimental  and  statistical  analysis 

3.1.  Single  cell  cold  start  experiments 


Table  1 

Varied  operating  parameters  for  start  up  experiments  at  — 10  °C 


Operating  parameter 

Values 

Applied  voltage  for  cold  start  up,  V  (mV) 

200/400/600 

Air  flow  rate,  ua ir  (ml  min-1) 

147-1000 

Hydrogen  flow  rate,  «h 9  (ml  min-1) 

20-158 

Cell  impedance  after  nitrogen  purging,  Zpurg ing  (m£2) 

70-1500 

Cell  impedance  after  applying  gases  directly  before 
cold  start,  Zstar t  (m£2) 

7-50 

Power  density  at  450  mV,  30  °C,  dry  gases, 
stoichiometry  1.5/3. 0, /7450  (mWcm-2) 

210-350 

After  this  reference  operation  the  gases  were  switched  to  dry 
nitrogen  (N2)  and  both  sides  of  the  fuel  cell  were  purged  with 
1 .0 1  min- 1  at  30  °C.  The  purging  process  continued  until  the  cell 
impedance  at  1000  Hz  increases  from  2.5  to  4.5  mQ  at  normal 
operation  conditions  to  values  larger  than  500  m£L  Then  N2  flow 
is  switched  off  and  the  cell  is  cooled  down. 

After  drying  and  freezing  of  the  cell,  the  start  up  experi¬ 
ments  were  carried  out  potentiostatically.  Gas  flow  rates  were 
kept  constant  as  well  as  the  cell  temperature.  When  the  graphite 
plates  reach  the  required  temperature,  dry  hydrogen  and  dry  air 
is  switched  on.  After  open  circuit  voltage  has  arisen  potentio- 
static  load  is  applied.  The  cumulated  charge  transfer  density 
S^cumW  serves  for  evaluation  of  the  cold  start  experiments  and 
is  calculated  by  (1) 

Sq,  cum(0  —  [  tdt  (1) 

Jo 

Sq, Cum(0  is  chosen  because  it  is  very  useful  to  compare  different 
cold  start  experiments  with  each  other.  The  higher  Sq?cum(t)  the 
better  is  the  cold  start  behaviour.  On  the  other  hand  Sq^cum(t) 
directly  correlates  with  the  amount  of  product  water  mH2o,prod 
given  by  Eq.  (2) 


^H20,prod 


}q,  cum 


(0 


2  F 


Mh2q 


In  (2)  mH2o,prod  is  the  area  specific  product  water  weight 
in  gem-2,  Sq,cum(t)  the  cumulated  charge  transfer  density  in 
Ccm-2,  F  the  Faraday  constant  96,485  C  mol-1  and  Mh2o  is 
the  molar  weight  of  water  in  g  mol-1 . 


3.2.  Experimental  cold  start  influence  parameters 

To  find  out  the  main  influencing  operating  parameters  15 
experiments  were  carried  out  at  constant  temperature  of  — 10  °C. 
In  Table  1  all  varied  operating  parameters  are  listed  and 
explained. 

In  a  second  experimental  series  the  influence  of  cell  compo¬ 
nents  was  studied.  In  Table  2  the  varied  fuel  cell  components 


Before  freezing  the  fuel  cell  a  reference  performance  mea¬ 
surement  at  0.45  V  with  dry  hydrogen  (stoichiometry  1.5)  and 
dry  air  (stoichiometry  3.0)  at  30  °C  was  recorded.  This  refer¬ 
ence  power  density  is  called  ^450.  Each  point  was  measured  at 
least  20  min,  gas  flows  were  adjusted  automatically  to  the  actual 
current. 


Table  2 

Varied  component  parameters  for 

— 10  °C  start  up  experiments 

Parameter 

Values 

MEA,  membrane  thickness 

18/25/35 

Gas  diffusion  layer 

GDL/MPL/OAC/OCC 

406 
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Fig.  1.  Open  channel  gas  diffusion  layer.  The  reaction  gas  can  diffuse  under  the 
landing  of  the  double  serpentine  flowfield  while  the  electrode  under  the  channels 
is  directly  exposed  to  the  gas  stream. 

are  listed.  Catalyst  coated  membranes  are  used  with  a  platinum 
catalyst  loading  of  0.4  mgpt  cm-2  on  both  electrodes.  As  gas 
diffusion  layer  (GDL)  carbon  cloth  was  used  with  a  thickness 
of  190  pim,  a  porosity  of  0.75  and  a  PTFE  content  of  10%. 
Experiments  with  standard  GDL  were  performed  as  well  as 
experiments  with  additional  micro  porous  carbon  layers  (MPL) 
on  both  electrodes. 

In  further  experiments  the  GDL  were  cut  in  a  way,  that  one 
electrode  surface,  either  the  anode  (OAC,  open  anode  channels) 
or  the  cathode  (OCC,  open  cathode  channels)  is  directly  exposed 
to  the  gas  stream  while  the  other  electrode  is  equipped  with  a 
standard  GDL.  Fig.  1  shows  an  open  channel  GDL. 

3.3.  Data  analysis  using  statistical  prediction  models 

The  statistical  software  Cornerstone®  (Brooks  Automation) 
fits  a  response  y,  e.g.  Sq,cum(t),  with  a  polynomial  term  consist¬ 
ing  of  a  constant  c,  the  values  jc*  of  the  p  influence  parameters, 
e.g.  applied  voltage  V,  multiplied  by  constant  factors  a^  In  (3) 
a  polynomial  term  is  shown,  where  linear  effects  are  taken  into 
account  in  the  first  term,  interaction  effects  between  each  param¬ 
eter  in  the  second  term  and  quadratic  effects  in  the  third  term 

p  p  p 

y  =  c  +  yfaixi  +  ai jxi j  +  'ffa»x2i  (3) 

i= 1  i,j  =  1  i=  1 

i  /  j 

The  fit  quality  of  the  regression  polynomial  is  evaluated  by  the 
stability  index  R 2  and  the  adjusted  stability  index  R 2d- .  R2  is  the 
relation  between  the  statistical  spread  of  the  model  and  the  statis¬ 
tical  spread  of  the  real  measured  values.  For  test  series  with  a  low 
number  of  experiments  R2d ■  is  better  suited,  because  it  accounts 
additionally  for  the  number  experiments  and  the  number  of  influ¬ 
ence  parameters.  The  values  of  R2  and  R2d ■  are  between  0  and 

1 .  For  serious  statements  R2d ■  should  normally  be  higher  than 
0.8  [5,6]. 

3.4.  Portable  fuel  cell  stack  cold  start  experiments 

To  investigate  the  starting  behaviour  of  fuel  cell  stacks 
two  test  series  were  performed  at  —10  and  —  20  °C.  Fuel  cell 


stacks  were  dried  by  N2  purging  whereby  the  stack  impedance 
increased  from  its  normal  value  of  20  to  100  mQ  and  above.  The 
gas  flow  rates  for  start  up  were  kept  constant  as  long  as  possi¬ 
ble  on  220  ml  min-1  H2  and  600  ml  min-1  air.  These  flow  rates 
allow  a  maximum  stack  current  of  more  than  5  A.  Measurable 
is  the  deviation  of  the  single  cell  voltages  from  the  average  cell 
voltage  and  the  time  until  the  fuel  cell  stack  temperature  reaches 
0°C. 

4.  Single  fuel  cell  results  and  discussion 

4.1.  Single  fuel  cell  drying 

Before  the  single  fuel  cell  is  cooled  down  below  0  °C  the  flow- 
field  must  be  free  of  water  droplets.  Therefore  drying  of  the  fuel 
cell  was  carried  out  until  the  cell  impedance  reaches  more  than 
500  m^.  If  the  purging  process  is  stopped  at  lower  impedance 
values  ice  blockages  in  the  gas  channels  were  observed.  This 
occurred  in  the  anode  and  also  the  cathode  flowfield,  detected 
by  increasing  inlet  pressures  when  hydrogen  and  air  are  applied. 
Because  of  the  fast  water  uptake  of  polymer  electrolyte  mem¬ 
branes  [7]  rest  humidity  from  the  electrodes  and  the  GDL  dif¬ 
fuses  back  into  the  membrane  when  the  drying  process  is  stopped 
and  impedance  could  drop  down. 

4.2.  Single  fuel  cell  operation  at  —10°  C  and 
parameterisation  of  results 

After  the  fuel  cell  is  dried  out,  cooled  down  to  — 10  °C  and 
Voc  has  arisen,  potentiostatic  load  is  applied.  The  resulting  cur¬ 
rent  is  decreasing  over  time  always  with  the  same  principle 
behaviour.  Fig.  4  shows  a  typical  current  decay  curve  over  time. 
The  fuel  cell  is  reacting  with  a  fast  increase  of  current  which 
reaches  within  the  first  2  min  its  maximum  and  then  decays. 
Also  illustrated  in  Fig.  4  is  the  cumulated  charge  transfer  den¬ 
sity  SqiCum(t).  To  predict  the  maximum  cumulated  charge  transfer 
density  Sq, max  a  mathematical  fitting  function  was  found  which 
parameterises  5^Cum(0-  This  exponential  function  is  calculated 
by  (4) 

=  \max(l  -  e-(f/*«)c)  (4) 

*63  is  the  time  when  Sq^ft)  reaches  63%  of  Sq,  max  and  in  combi¬ 
nation  with  the  variable  c  it  describes  the  slope  of  the  curve.  This 
mathematical  fitting  curve  is  also  shown  in  Fig.  2.  The  quality 
of  a  cold  start  experiment  is  evaluated  by  Sq> max  as  response 
parameter  in  the  statistical  regression. 

Because  the  pressure  drop  of  hydrogen  and  air  over  the  flow- 
field  remains  stable,  we  assume  that  most  of  the  product  water 
is  not  transported  into  the  channels,  but  freezes  in  the  porous 
structure  of  the  GDL  and  the  electrode.  As  result  the  reaction  is 
inhibited. 

To  analyse  the  initial  starting  behaviour  the  current  decay 
is  referred  to  the  cumulated  charge  transfer  density  Sq,cum(t). 
Fig.  3  illustrates  current  density  i  versus  Sq^cum(t)  for  the  same 
experiment  as  shown  in  Fig.  2.  For  statistical  evaluation  it  is 
necessary  to  parameterise  these  current  density  curves.  A  simple 
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Fig.  2.  Decay  of  current  density  and  cumulated  charge  transfer  over  time  at  a 
cell  voltage  of  200 mV,  cell  temperature:  — 10  °C,  gas  flows:  600mlair  min-1, 
120  mlH2  min-1 ,  cell  impedance  at  start  1 1 .2  m£2  and  a  reference  power  density 
P450  at  30  °C  of  250 mW cm-2. 


way  is  to  fit  the  almost  straight  current  curve  sections  by  straight 
lines  as  shown  in  Fig.  3.  This  fit  is  created  by  manual  selection 
of  five  characteristic  points,  which  describe  the  principle  curve 
behaviour. 


4.3.  Statistical  analysis  of  influence  of  operational 
parameters  on  Sq>max 

A  linear  regression  analyses  was  performed  to  predict 
the  main  influencing  operational  parameters  on  Sq,max.  Three 
parameters  were  found  which  have  a  significant  impact  on  the 
value  of  S^max-  The  largest  effect  has  the  reference  power  den¬ 
sity  P450.  The  influence  of  the  air  flow  rate  u air  and  the  membrane 
humidity  indicated  by  Zstar t  is  lower.  The  coefficient  correlation 
between  these  three  parameters  is  lower  than  0.4,  this  means 
that  they  are  independent.  Changes  in  the  reference  cell  power 
density  ^450  are  not  due  to  potential  changes  of  the  membrane 
water  uptake.  During  reference  operation  no  significant  changes 
of  the  cell  impedance  at  1000  Hz  were  observed. 

Fig.  4  illustrates  the  regression  result  for  Sq, max  as  predicted 
response  graph.  On  the  y-axis  the  predicted  response  from  the 
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Fig.  3.  Current  density  decay  over  cumulated  charge  transfer  density  Sq,cum(t) 
and  manual  fitted  five-point  curve.  Cell  voltage:  200  mV,  cell  temperature: 
— 10  °C,  gas  flows:  600  mlair  min-1, 120mlH2  min-1,  start  impedance:  11.2m£2, 
reference  power  density:  250  mW  cm-2. 
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Fig.  4.  Predicted  S^max  as  function  of  the  dominating  three  operational  parame¬ 
ters.  Linear  regression  of  all  operational  parameters  with  a  regression  goodness 
of  R^d .  =  0.87  and  a  root  mean  squared  error  of  4.3  C  cm-2. 


regression  polynomial  is  shown  as  function  of  the  values  of  J7450 
(250mWcm-2),  u air  (600  ml  min-1)  and  Zstar t  (40  mQ).  The 
predicted  value  for  S^max  is  47.99  C  cm-2  with  a  standard  devi¬ 
ation  of  ±12.9  C  cm-2. 

This  result  shows,  that  two  operational  factors  can  be  adjusted 
by  the  operator  to  achieve  a  high  load  transfer  rate.  The  mem¬ 
brane  should  be  dry,  indicated  by  a  higher  Zstar t  and  the  air  flow 
rate  should  be  also  as  high  as  possible.  Both  can  be  explained 
by  the  influence  product  water.  The  dryer  the  membrane,  the 
higher  is  its  water  uptake  and  the  higher  the  air  flow,  the  bet¬ 
ter  is  the  removal  of  product  water  on  the  cathode  side.  The 
third  and  most  influencing  parameter  is  J7450.  This  reference 
power  density  is  strongly  depending  of  the  degradation  status 
of  the  cathode  electrode.  When  the  electrode  performance  is 
low  because  of  degradation,  ^450  is  also  lower.  This  means  that 
during  the  potentio static  cold  start  the  output  power  and  there¬ 
fore  the  heat  production  is  reduced  too.  This  is  why  in  the  case 
of  low  P450  the  electrode  temperature  is  lower  compared  to  cases 
with  a  higher  P450.  This  could  lead  to  a  faster  flooding  and  ice 
formation  in  the  electrodes. 


4.4.  Statistical  analysis  of  the  influence  of  components  on 

Sq,  max 


The  influence  of  components  is  strongly  dominated  by  the 
GDL.  The  experiment  space  for  the  GDL  porosity  is  too  huge 
and  therefore  a  statistical  regression  analysis  is  not  able  to  show 
significantly  the  influence  of  the  also  investigated  membrane 
thickness.  Sq,mSLX  is  mainly  influenced  by  the  GDL,  ^450  and  the 
air  flow  u air.  The  goodness  of  the  regression  analysis  is  not  as 
high  as  for  the  operational  parameters  with  7?2d ■  =0.56  and  root 

mean  squared  error  of  17.7  C  cm-2. 

Fig.  5  illustrates  the  values  of  S^max  grouped  by  the  GDL 
configuration.  A  more  open  structure  of  the  GDL  leads  to  higher 

Sq,max  • 

An  interesting  fact  is  that  an  open  anode  GDL  structure  led 
to  the  highest  values  for  S^max-  This  indicates  that  not  only 
the  cathode  electrode  is  responsible  for  the  inhibit  of  the  reac¬ 
tion.  Even  at  low  temperatures  a  strong  water  permeation  to 
the  anode  side  seems  to  exist,  otherwise  the  impact  of  the 
open  anode  channel  structure  could  not  be  explained.  Further 
investigations  should  be  performed  to  find  out,  if  primarily 
the  cathode  electrode  freezes,  as  proposed  by  Kagami  et  al. 
[7],  or  if  the  anode  electrode  could  also  be  affected  by  ice 
formation. 
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Table  3 


Coefficients  ai  for  polynomial  prediction  model 


(1)  /start 

(2)  5g5peak;jStart 

(2)  /peak 

(3)  ‘S'^peal^end 

(4)  Sq, linear 

(4)  ^linear 

(5)  5)^max 

(mAcm-2) 

(niA  cm-2) 

(mA  cm-2) 

(mAcm-2) 

(C  cm-2) 

(mA  cm-2) 

(mA  cm-2) 

Constant 

-977.8981 

-22.4874 

144.9494 

7.4862 

3.9557 

-1.9906 

-26.0127 

V 

2.9383 

-0.2209 

-0.0180 

-0.0056 

^air 

0.9819 

0.0229 

-0.4419 

0.0104 

-0.0024 

0.0582 

0.0247 

WH2 

0.6911 

-0.6934 

Ztart 

-36.3194 

0.3027 

-3.1301 

0.3954 

0.4184 

0.4772 

P450 

3.8487 

0.1402 

0.1094 

0.1663 

Wair^H2 

Vfiair 

-0.0014 

-0.0008 

-0.0001 

TZstart 

0.0158 

0.0036 

P450V 

Zstart 

-0.0104 

-0.0078 

0.0107 

^ai  restart 

-0.0123 

0.0004 

0.0004 

P450  wair 
P45()Zstart 

0.0888 

0.0021 

0.0009 

Fig.  5.  S^max  as  function  of  GDL  configuration. 


4.5.  Influence  of  operational  parameters  on  current  decay 
and  interpretation 


Fig.  6.  Statistical  based  fit  (filled  symbols:  ♦  A)  of  four  current  decay 

curves  at  —10  °C  at  different  operating  conditions  as  described  in  the  legend. 


To  find  out  which  operational  parameters  affect  the  shape 
of  the  current  decay  curve  linear  regression  analysis  was  made 
for  each  point  (see  Fig.  6).  Interactions  of  the  main  parameters 
(£/,  wair,  uw2,  Zstart,  Paso)  are  included.  The  goodness  of  the 
regression  analysis  is  indicated  by  R ^  which  is  for  all  points 
between  0.855  and  0.998,  except  for  point  5  (S^max),  where 
R^d j  is  0.640.  Table  3  shows  all  coefficients  at  for  the  statistic 
prediction  polynomial. 

As  an  example  the  regression  model  calculates  Sq, linear  (point 
4)  with  a  polynomial  expression,  given  by  the  constant  c  and 
factors  at  multiplied  by  the  adjusted  values  influence  parameter, 


in  this  example  ua jr  and  Zsta rt- 
sq, linear  =  3.9557  -  0.0024wair  +  0.4184Zstart 
T0.0004wair  Zstart 

The  quality  of  the  model  is  demonstrated  in  Fig.  6.  Points  1-4 
are  well  fitted,  point  5  not  because  of  the  low  Rfd y . 

The  statistic  based  model  estimates  the  dominating  factors 
for  each  curve  point.  In  Table  4  the  influence  is  shown. 

The  regression  shows  that  the  initial  start  current  z'start  is 
mainly  depending  on  the  membrane  humidity  Zstar t  and  the 


Table  4 


Estimated  influence  of  main  parameters  on  current  decay  curves 


Parameter 

(1)  4 tart  (%) 

(2)  5g5peak, start  (%) 

(2)  ipeak  (%) 

(3)  ^pea^end  (%) 

(4)  5"^^ linear, start  (%) 

(4)  ^linear, start  (%) 

(5)  \max  (%) 

V 

34 

— 

40 

25 

— 

59 

— 

^air 

19 

68 

92 

5 

13 

20 

5 

WH2 

— 

69 

— 

4 

— 

— 

— 

Zstart 

43 

36 

32 

87 

85 

— 

26 

P450 

12 

— 

34 

— 

— 

10 

26 

M.  Oszcipok  et  al.  /  Journal  of  Power  Sources  154  (2006)  404-411 


409 


operation  voltage  V.  We  assume  that  humidification  of  the 
membrane  and  the  anode  electrode  sets  in  by  product  water  until 
the  current  density  reaches  its  peak  plateau  /peak-  The  mem¬ 
brane  reaches  maximum  protonic  conductivity  at  SVy, peak, end- 
Now  product  water  starts  to  fill  mainly  the  pores  of  the  elec¬ 
trodes  and  the  current  density  drops  down  very  fast.  The  level  to 
which  it  drops  down  /linear, start  is  strongly  dominated  by  the  oper¬ 
ating  voltage  V.  Lower  voltage  is  connected  with  higher  waste 
heat  production  at  the  electrode  resulting  in  a  higher  /linear, start 
because  more  product  water  stays  in  a  liquid  form  and  can  be 
carried  out  by  the  air  flow. 

This  interpretation  is  supported  by  the  fact,  that  all  charge 
transfer  related  statements  are  highly  influenced  by  the  drying 
status  of  the  cell  whereas  all  current  related  statements  can  be 
interpreted  through  heat  production  and  are  significantly  influ¬ 
enced  by  the  cell  performance  ^450,  the  cell  voltage  V  and  the 
air  flow  u^. 

These  results  show,  that  the  presented  current  decay  curves 
are  very  complex  and  depend  on  a  variety  of  influence  param¬ 
eters.  Existing  physical  fuel  cell  models  cannot  directly  be 
applied.  A  combination  of  several  model  approaches  is  seen 
to  be  promising  [7,8]. 


4.7.  Conclusions  of  isothermal  single  cell  experiments 

Our  statistical  based  model  estimates  the  influencing  parame¬ 
ters  of  isothermal  cold  start  experiments  of  single  fuel  cells.  The 
results  support  the  work  of  Kagami  et  al.,  that  a  high  air  stream 
delays  the  reaction  decay.  Additionally  the  results  show,  that  the 
influence  of  the  starting  impedance  and  the  overall  cell  perfor¬ 
mance  at  30  °C  have  an  important  influence  too.  We  assume  that 
a  product  water  flux  from  the  cathode  to  the  anode  influences 
the  initial  cold  start  behaviour  in  particular  and  also  the  time 
until  the  cell  reaction  is  disrupted.  It  is  found  that  the  overall 
fuel  cell  performance  at  450  mV  and  30  °C  is  strongly  affected 
by  the  cold  start  experiments  and  has  also  a  strong  influence  on 
the  cold  start  ability.  The  influence  of  the  heat  production  con¬ 
nected  with  applied  cell  voltage  for  cold  start  is  low  at  isothermal 
operation.  Investigation  of  the  component  influence  showed,  that 
a  high  GDL  porosity  could  enable  a  better  water  removal  and 
therefore  a  higher  cumulated  charge  transfer. 

5.  Results  and  discussion  of  fuel  cell  stack  experiments 

5. 1.  Cold  start  experiments  at  —10°  C 


4.6.  Performance  loss  by  cold  start  experiments 

In  Fig.  7  the  reference  measurements  at  450  mV  before  each 
cold  start  experiment  is  shown.  A  performance  loss  is  obvi¬ 
ous  for  each  cell  assembly.  It  is  assumed  that  reduction  of 
the  effective  catalyst  surface  could  be  made  responsible  for 
this  decay.  Oh  et  al.  reported  from  catalyst  layer  degrada¬ 
tion  because  of  volume  expansion  during  the  phase  change 
of  liquid  water  to  ice.  The  effective  active  surface  is  reduced 
and  the  pore  size  distribution  in  the  catalyst  layer  changes 
[9] .  Additionally  slight  changes  of  membrane  properties  could 
also  affect  the  performance  loss.  McDonald  et  al.  found  out 
that  membrane  swelling  by  water  uptake  is  reduced  during 
freeze  thaw  cycling,  but  the  electrochemical  properties  were  not 
influenced  [10]. 


In  a  first  test  series  the  climate  chamber  was  adjusted  to  a  tem¬ 
perature  of  — 10  °C.  Three  start  up  experiments  were  performed 
were  the  starting  impedances  were  between  104  and  660  m^. 
The  only  regulation  parameter  was  the  applied  load  to  the  fuel 
cell  stack. 

Fig.  8  shows  the  first  start  up  experiment  at  1  A.  In  the  graph 
the  average  single  cell  voltage  of  the  stack  with  the  minimum 
and  maximum  single  cell  voltages  can  be  seen  as  well  as  stack 
temperature  and  stack  current.  The  waste  heat  <2waste  is  also 
shown.  It  is  calculated  by  (5)  using  the  enthalpy  of  formation 
per  mole  —  A/if,  the  Faraday  constant  F ,  the  average  cell  voltage 
Vceii,  the  number  of  cells  n  and  the  stack  current  /stack,  assuming 
that  the  product  water  is  condensed  back  to  liquid  [11] 


Q Waste 


/  —  A/if 

v  2F 


n  /stack 


Fig.  7.  Performance  loss  at  450  mV  during  four  different  cold  start  series.  Ref¬ 
erence  measurements  were  performed  at  30  °C  cell  temperature,  dry  reaction 
gases  and  constant  stoichiometries  of  1.5  for  H2  and  3.0  for  air.  The  minimum 
measuring  time  for  each  data  point  was  20  min. 
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Fig.  8.  Start  up  experiment  of  a  six-cell  fuel  stack  at  — 10  °C.  The  impedance  at 
start  was  660  m£2,  gas  flow  rates  were  220  ml  min-1  H2  and  600  ml  min-1  air. 
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Fig.  9.  Start  up  experiment  of  a  six-cell  fuel  stack  at  — 10  °C.  The  impedance  at 
start  was  230  m£2,  gas  flow  rates  were  220  ml  min-1  H2  and  600  ml  min-1  air. 


Fig.  11.  Start  up  experiment  of  a  six-cell  fuel  stack  at  —20  °C.  The  impedance 
at  start  was  1 20  mQ,  gas  flow  rates  were  220  ml  min- 1  H2  and  600  ml  min- 1  air. 


In  the  beginning,  the  single  cell  voltages  spread  very  strong 
between  260  and  800  mV.  In  the  following  the  voltages  converge. 
Responsible  is  the  very  dry  state  of  the  stack  (Zsta rt  =  660  m£2) 
which  is  later  humidified  by  product  water.  The  stack  tempera¬ 
ture  raised  up  to  —5  °C  within  the  first  7  min.  Then  the  load  was 
increased  to  1.5  and  2  A.  After  17  min  the  stack  temperature 
reached  0  °C. 

In  Fig.  9  the  start  up  is  shown  for  a  current  increase  to  2  A. 
The  start  impedance  was  230  mQ  and  therefore  the  single  cell 
voltages  are  distributed  more  equally.  The  stack  temperature 
increased  in  6.5  min  to  0  °C. 

Fig.  10  illustrates  an  experiment  with  a  current  increase  up 
to  2.8  A.  The  starting  impedance  was  100  mQ.  A  temperature  of 
0  °C  was  reached  within  4  min. 


5.3.  Discussion  of  cold  start  experiments  on  portable  fuel 
cell  stacks 

At  temperatures  of  — 10  °C  and  higher  it  was  possible  to  start 
a  portable  fuel  cell  stack  passive.  In  all  experiments  the  stack 
temperature  increased  to  operating  temperatures  over  40  °C  at  a 
constant  ambient  temperature  of  — 10  °C.  The  impedance  seems 
to  have  also  in  fuel  cell  stacks  a  strong  influence  on  the  start  up 
behaviour.  If  the  stack  is  too  dry,  the  behaviour  of  the  single  cells 
differs  strongly  and  does  not  allow  an  increase  of  load  which 
would  lead  to  a  higher  waste  heat  production  of  the  stack.  At 
— 20  °C  the  single  cells  react  very  sensitive  when  load  is  applied. 
With  the  described  experimental  set-up  a  successful  cold  start 
up  was  not  possible. 


5.2.  Cold  start  experiments  at  —20°C 

In  a  further  test  series  the  temperature  was  lowered  to  —  20  °C. 
The  gas  flow  rates  were  kept  constant  and  the  only  start  up  regu¬ 
lation  was  done  by  adjustment  of  the  applied  load.  Fig.  1 1  shows 
a  cold  start  experiment  at  —  20  °C.  The  cell  reacts  more  sensitive, 
when  load  is  applied  than  at  — 10  °C.  It  was  not  possible  to  reach 
0  °C  with  passive  methods  whether  in  this  case  or  other  —  20  °C 
cold  start  experiments. 
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Fig.  10.  Start  up  experiment  of  a  six-cell  fuel  stack  at  —10  °C.  The  impedance 
at  start  was  100  mQ,  gas  flow  rates  were  220  ml  min- 1  H2  and  600  ml  min-1  air. 


6.  Conclusions 

Substantial  experiments  of  the  cold  start  behaviour  of  fuel 
cells  are  presented,  beginning  with  basic  investigations  on  single 
cells  and  moreover  cold  start  up  experiments  of  portable  fuel  cell 
stacks  in  a  climate  chamber.  An  empirical  based  statistical  model 
was  introduced  which  describes  the  initial  cold  start  behaviour 
of  single  cells  in  good  correlation  to  experimental  results.  The 
main  influencing  operational  parameters  were  identified.  Results 
of  the  single  cell  experiments  could  be  successfully  adapted  to 
start  up  portable  stack  from  — 10  °C.  The  load  level  must  be  high 
for  a  high  waste  heat  level,  the  air  flow  rate  should  be  high  to 
blow  out  product  water  and  for  successful  cold  start  up  in  fuel  cell 
stacks  the  membrane  should  not  be  too  dry.  Degradation  effects 
in  single  cells  were  observed  resulting  in  lower  performance  at 
30  °C. 
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